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By mapping tlie relativistic version of the Chern-Simons-Landau-Ginzburg theory in 2+1 dimen- 
sions to the 3D lattice Villain x-y model coupled with the Chern-Simons gauge field, we investigate 
phase transitions of Chern-Simons bosons in the limit of strong coupling. We construct algebraically 
exact duality and flux attachment transformations of the lattice theories, corresponding to analogous 
transformations in the continuum limit. These transformations are used to convert the model with 
arbitrary fractional Chern-Simons coefficient a to a model with a either zero or one. Depending 
on this final value of a, the phase transition in the original model is either in the universality class 
of 3D x-y model or a "fermionic" universality class, unless the irrelevant corrections of cubic and 
higher power in momenta render the transition of the first order. 



INTRODUCTION 



As one changes the external magnetic field, the Hall 
conductance of a system of two-dimensional electrons 
jumps from one quantized value to another, while the lon- 
gitudinal conductance displays a peak. The zero temper- 
ature localization length diverges in transition points, in- 
dicating the second order quantum phase transitions. Re- 
markably, experiments ||l|-^ show that these phase tran- 
sitions are universal, i.e. the critical exponents governing 
the divergence of the localization length are the same for 
integer and fractional transitions. It has also been found 
Q that the critical resistances for integer and fractional 
phase transitions are universal as well. Finally, a bril- 
liant recent experiment by Shahar et al demonstrates 
precise duality mapping of non-linear current-voltage 
characteristics between the insulating and the fractional 
quantum Hall sides of the transition, where carriers are 
respectively electrons and fractionally charged quasipar- 
ticles. One of the most challenging problems in the quan- 
tum Hall effect is to understand the universality of the 
transitions. 

Apparently, there is a conceptual difference between 
the physics of the integer and the fractional quantum 
Hall regimes. In the former case the quasiparticles are 
fermions, their localization transition is believed to be a 
one-particle problem. The physics is understood qua- 
siclassically as quantum tunneling-assisted percolation 
1^-^, while the non-linear sigma model with topological 
term 1^ provides the field-theoretical formalism to study 
this problem. On the other hand, the excitations above 
the FQHE ground state, given with extreme precision 
by inherently many-body Laughlin's wavefunction |10|, 



have the infinite-range statistical interaction that seem 
to invalidate the usual one-particle localization theory. 
Field-theoretically, this interaction can be described in 
terms of the Chern-Simons gauge field coupled to either 
bosons |l|,|l|, or fermions ||-|l§. 

Using the bosonic description, called the Chern- 
Simons-Landau-Ginzburg (CSLG) theory of the 
quantum Hall effect, Kivelson, Lee and Zhang (KLZ) pro- 
posed |l^ a global phase diagram of the quantum Hall 
effect. The original problem of interacting electrons in 
large magnetic field is mapped to the problem of bosons 
in zero or weak uniform magnetic field, interacting via 
the Chern-Simons gauge field. The information about 
the filling factor is contained only in the coefficient of the 
Chern-Simons term. The electromagnetic response func- 
tions in the CSLG theory can be expressed solely in terms 
of response functions of the bosonic field and, therefore, 
the different quantum Hall transitions are in the same 
universality class if the critical exponents of the boson 
superfluid to insulator phase transition are independent 
of the Chern-Simons coefficient. 

On the basis of the CSLG theory in the RPA approx- 
imation KLZ argued that the critical exponents of Hall 
transitions are universal; they predicted values for critical 
conductances that have been recently confirmed experi- 
mentally [Q. However, since the Chern-Simons term is a 
marginal operator by naive power counting, it could in 
principle change the critical exponents. If this happens, 
the transition will no longer be universal. Therefore, it is 
highly desirable to check the idea of universality in some 
specific field-theoretical calculation. Unfortunately, it is 
very hard to treat the full interacting problem in the 
presence of disorder. For this reason the analysis is of- 
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ten restricted to toy models that nevertheless share some 
essential features with real quantum Hall systems. 

The most spectacular feature of the quantum Hall sys- 
tem is the symmetry of the phase diagram, and there 
are several field-theoretical models with similar symme- 
try. A broad class of such models has been introduced by 
Shapere and Wilczek jlj] who generalized the construc- 
tion of Cardy and Rabinovici |19 2^ and constructed 4 
and 2 dimensional self-dual models consisting of mutually 
dual discrete Abelian gauge lattice models coupled via 
the theta term. The symmetry of similarly constructed 
2-1-1 dimensional model with Chern-Simons term was 
analyzed by Rey and Zee [Q. Liitken and Ross [ p2|j23| 
discussed the symmetry of fixed points of self-dual 2- 
dimensional clock model and argued that the correspond- 
ing group SL{2,Z) describes the symmetry of all fixed 
points of the quantum Hall system. This statement is 
analogous to the Law of Corresponding states im- 
fortunately, so far no derivation of the field-theoretical 
realization of this symmetry group is known. 

The question whether phase transitions in the system 
of Chern-Simons bosons can even in principle be univer- 
sal is so general that the answer for virtually any such 
model would be interesting. Particularly, one can study 
the phase transition driven by the amplitude of external 
periodic potential at commensurate filling instead of that 
driven by disorder. It is known that the system of 
bosons at fillings commensurate with external periodic 
potential can be mapped to the x-y model. In the quan- 
tum Hall problem one arrives p5|] at the x-y model mini- 
mally coupled with the Chern-Simons term in zero aver- 
age magnetic field. Within this model one can naturally 
address the question whether the order-disorder transi- 
tion is altered by coupling with the Chern-Simons term 
and, if so, whether there is a new universality class, or 
the phase transition continuously depends on the Chern- 
Simons coupling. Even though this is an artificial model, 
its successful analysis will undoubtedly shed some light 
on the universality of the real quantum Hall phase tran- 
sitions. 

Wen and Wu |25| used the representation of the the x-y 
model in terms of the relativistic complex-valued scalar 
field with quartic interaction, especially convenient for 
perturbational expansion. They also used the 1/7V ex- 
pansion, known to work for ordinary bosons, to regularize 
the perturbational series in three dimensions. Within the 
second loop approximation Wen and Wu found that al- 
though the phase transition remains of the second order, 
the critical exponents of the x-y model are continuously 
changed by the Chern-Simons coupling. 

Soon after that, we pointed out that the 1/N ex- 
pansion artificially suppresses the gauge fiuctuations in 
this model, and the result of Wen and Wu is not likely to 
hold in the physical limit of iV = 1 . Our alternative cal- 
culation was performed in the vicinity of nominal tri- 
critical point of the theory, where the quartic interaction 
is absent and one has to keep the six-field coupling. The 



effective theory in the vicinity of the tricritical point is 
renormalizable in three dimensions and allows systematic 
loop expansion, exactly as the theory with four-field in- 
teraction is renormalizable in four dimensions where the 
quartic coupling is dimensionless. The presence of the 
Chern-Simons coupling, also dimensionless in 3D, does 
not introduce any classically divergent quantities, but it 
has a dramatic effect already at the second loop level, 
changing the transition from the second to the first order. 
This is similar to the statement |^ that the phase tran- 
sition in clean type I superconductors becomes weakly 
first order due to the fiuctuations of the electromagnetic 
field. In our case, since the theory without the quartic 
interaction is massive, it remains massive in some finite 
region of values of the quartic coupling. Thus, within the 
perturbative region of the relativistic complex scalar field 
model, the transition is of the first order even for some 
non-zero values of the quartic coupling. 

This conclusion obtained in the vicinity of the Gaus- 
sian fixed point of the scalar sector was not in formal 
contradiction with the result of Wen and Wu obtained 
near the strongly coupled fixed point (formally accessed 
with the 1/A'' expansion,) but it reopened the question 
about the universality of phase transitions of scalar fields 
coupled with the Chern-Simons field. The results of both 
calculations indicate that the Chern-Simons interaction 
is a relevant perturbation to the scalar theory with strong 
short-range repulsion and some kind of non-perturbative 
analysis is necessary. 

The goal of this paper is to address the question of uni- 
versality of phase transitions in relativistic Chern-Simons 
bosons without the limitations of the perturbation the- 
ory. We access the strong coupling limit of the Chern- 
Simons bosons using the Villain form of the x-y model 
minimally coupled with the Chern-Simons gauge field. 
For this model we formulate the exact duality and fiux 
attachment transformations with usual properties in the 
continuum limit and build the Haldane-Halperin ^8 29 
hierarchy of models related by these non-local transfor- 
mations. We show that in the absence of external mag- 
netic field the universality class of phase transition in 
the strong-coupling limit is determined by the Chern- 
Simons coupling a = p/q. When either p or g is even, 
the phase transition is in the universality class of the 
order-disorder transition in the usual three dimensional 
x-y model. When both p and q are odd, the system in the 
vicinity of critical point (if any) is equivalent to x-y model 
with the Chern-Simons coupling a = 1, corresponding to 
the Fermi statistics. We believe that our result is exact 
and apply at least for fractions with small enough de- 
nominators q. With increasing denominators the terms 
formally irrelevant near the x-y critical point grow, ef- 
fectively limiting the number of fractions with expected 
universality. 

Unlike the previous works [|30|-|35t analyzing lattice 
Chern-Simons models, we derive algebraically exact fiux 
attachment and duality transformations. Thus we avoid 
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the problem of relevance or irrelevance of the higher or- 
der in momenta terms appearing in different lattice def- 
initions of the Chern-Simons coupling. Instead, we con- 
struct the exact transformation to the theory with zero 
a, equivalent to the usual 3D x-y model with finite lattice 
temperature. For this theory the algebra of local critical 
operators is known, and the formal irrelevance of extra 
terms does not require additional proof. 

The paper is organized as follows. We derive the ex- 
act duality transformation for the Villain form of the x-y 
model in the Section ^ and the exact flux attachment 
transformation in the Section These transformations 
have correct nai've continuum limit and work for lattice 
theories with gauge coupling of most general form. In the 
third section we construct the sequence of such transfor- 
mations leading to Villain x-y model with the Chern- 
Simons coupling either zero or one and utilize the known 
properties of the usual x-y model to build the phase di- 
agram of the original model. 



I. GENERALIZED DUALITY 

Particles in superfluid repel each other at short dis- 
tances. At very small temperatures most of the particles 
are in the same quantum state, characterized by the con- 
densate wavefunction = ^expiO, where both the 
amplitude tp and the phase are real-valued functions. 
Slow phase rotations have extremely small energy, imply- 
ing the existence of a linear sound mode. 

Another important kind of excitations in superfluid are 
vortices in two dimensions, or vortex lines in three dimen- 
sions. They are characterized by non-trivial phase of the 
condensate wavefunction gained along any surrounding 
contour. Far enough from the center of the vortex the 
gradient of phase is small and the superfluid density tp"^ 
is close to its non-perturbed value. Although the reduced 
density in the core region near the center of the vortex 
costs some energy, the total energy of a single vortex 
is mostly determined by the twisted phase in the area 
outside the core and any two vortices have a long range 
logarithmic interaction. 

Quite differently, vortices in superconductors have only 
short-range interactions (just like particles in superfluid): 
the extra phase is easily screened by the vector poten- 
tial. Instead, the conserved particle currents in nor- 
mal phase interact logarithmically through their mag- 
netic fields just like vortices in superfluid. This analogy 
can be formulated more precisely as the duality between 
the type II superconductors in London limit where the 
penetration length of the magnetic field is large com- 
pared to the coherence length and the strong-coupling 
limit of the Ginzburg-Landau model. In this limit of 
strong coupling one can completely neglect the density 
fluctuations, considering only the phases On as the only 
degrees of freedom. The core energy of the vortices can 
be regularized by redefining the theory at the lattice. 



The resulting x-y model is a collection of classical spins 
s„ = (cos 0„, sin 0„) with nearest-neighbor interaction of 
the form s^Sj = cos(6'i — 9j). It turns out that the re- 
verse transformation from the x-y model to the strongly 
coupled superfluid is also possible [|6|. This two-way 
analogy has been confirmed by comparing the results of 
numeric computations in the lattice x-y model with cor- 
responding analytical results for the cj)'^ theory. 
The partition function of the x-y model 



^x-y — / t; — 



\ - cos I 



2n ^"P^ 



-On) 



T 



(1) 



is periodic with respect to the phases On and, at least in 
the strong coupling limit where the lattice temperature 
p7| T is small, is mainly determined by the vicinity of 
maxima of the expression in the exponent. The shape of 
these maxima can be simulated using the Villain form of 
the x-y model 



exp ■ 



cos 



-On) 



T 



exp - 



,-0n-2TTmnf,] 

2T 



originally proposed by Berezinskii ||38|| and later indepen- 
dently and in greater details investigated by Villain . 

Peskin Q proved that the duality transforms the Vil- 
lain x-y model into the so-called frozen superconductor, 
or the zero temperature limit of the Villain x-y model 
coupled to the Maxwell field. This rigorous analysis has 
been extended by Kleinert |^ who analyzed the lattice 
superconductor in wide region of parameters using the 
Villain approximation and achieved quantitative agree- 
ment with numerical simulations. These model com- 
putations eventually lead [^^-Q to understanding of the 
complete phase diagram of superconductors in the dual 
representation. Later Lee and Fisher |45|j4^ ] and Lee 
and Zhang applied this duality transformation to the 
anyon superconductivity and the quantum Hall effect. 

We are interested in a very similar class of models, 
namely the strong coupling limit of relativistic Chern- 
Simons bosons. The idea is to follow the successful ex- 
ample of superconductor and examine non-perturbative 
properties of this model at the lattice using the Villain 
form of x-y model coupled to the Chern-Simons gauge 
field. There are several |^,^,^,^ different gauge- 
invariant definitions of the lattice Chern-Simons cou- 
pling, and, in order to keep the analysis as universal as 
possible, we consider the three-dimensional Villain x-y 
model 

/t^E exp- (^^^"-^:i;-^-"^)' , (2) 



na „ m— — oo 

" — TV 



2T 



Z = I exp (^-^A„M^::::.A.v) ZlAlUdAnp, 



(3) 



at the cubic lattice, minimally coupled to the gauge field 
An^ with some generic, possibly non-local kernel /C^ 
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or /C'^'^(k) = /C(k) in momentum representation. By the 
generic Chern-Simons coupling we shall imply that in 
the continuum limit K. must be an antisymmetric matrix 
linear in momenta. 

The summation over links /i in is performed 

only in the positive direction ft; the phases 6n are de- 
fined in the nodes, while both the integer-valued field 
ninfi and the gauge field are defined at the links 
of three-dimensional cubic lattice. We use the notations 
A^/„ = fn+p. - fn, A^/„ = fn- fn-fi for the forward 
and backward lattice differences respectively, and imply 
the summation over repeated indices. 

It will be convenient to specify the Fourier expansion 
of vector fields in the form 



1 



(2^ 



(4) 



The integration by 0„ can be done after regrouping the 
terms in the exponent. 



(8) 



the resulting constraint being merely the lattice version 
of the equation V • b = indicating that the field b is a 
pure curl. 



PP- 



(9) 



So far the only difference with the duality transfor- 
mation fioll for the pure x-y model is the presence of 
the product —ibn^Anfj, in the exponent of This extra 
term in the action is responsible for important additional 
symmetry. Indeed, the original gauge field in the full par- 
tition function (0) can be integrated away, 



(no summation in /i!) to emphasize their location at the 
links of the lattice. With this definition the gradient 
A^/„ = A^/„_|_^ or P^/(k) in momentum representation 

is a vector, while the divergence A^Anp, = A^An-jifi or 
P^j4p(k) is a scalar field, as one would expect from the 
geometrical interpretation of the lattice fields. We shall 
use the vector = 2sin(fc^/2) to define different func- 
tions of lattice momentum. For example, the kernel of 
the local in the gauge field Chern-Simons term used in 
refs. [||,|l]j3|] can be written as /C^^ = e'"'PPpQa/27:a 
with Qo = cos^^fc^/2, and in the Section || we shall 
define the Chern-Simons term with the kernel ICf^^ = 
sP-'^pPp/2TTaQo providing for the local coupling between 
the gauge-invariant scalar currents on the lattice. 

Let us treat the action Z[A] as the Villain x-y model 
in the presence of some external gauge field A, and follow 
the derivation of the duality transformation for this 
model. While the integration in phases 0„ is performed 
over restricted intervals, the gauge transformation 



27riV„ 



(5) 



does not change the integrand of (g). One can extend 
the integration region in On to the infinite interval by 
simultaneously constraining the field m„^ to avoid over 
counting. This can be done, for example, by applying 
the usual gauge condition 



rrinfj.) . 



(6) 



Now introduce an auxiliary field by writing 



(7) 



X exp 



T 

'—bl^ + ibn^ (Ap6'„ - - 27rTO„p) 



71 fl 

X exp 



27r 



1/2 



n 



2] J2^^^^-^^ (r<5,.+/c;J)6^(k) 



(10) 



and, since the field b is transverse, it is clear that the 
partition function Z depends only on the transverse com- 
bination 



(11) 



indicating the exact equivalence of any two models with 
the parameters satisfying the relationship 



(12) 



where = P^P^ denotes the transverse part of the 

Kronecker symbol and P^ = Pi//|P| is unit vector. This 
reparametrization changes both the gauge and scalar cou- 
plings; it is an exact symmetry of the partition function 
of the Villain x-y model or any correlators that do not 
involve the gauge field directly. These are the only physi- 
cal correlators as long as we treat the Chern-Simons field 
as an auxiliary field needed to define the fractional statis- 
tics for the particles. Since the lattice temperature 
T measures the strength of the local repulsion, the rela- 
tionship ( |l^ implies that a part of this repulsion can be 
mediated by the gauge field if its propagator has non-zero 
symmetric part. 

One may further verify that the system averages taken 
in the presence of any number of vortex-antivortex pairs 
expi(0„ — On'), as well as the monopoles dual to them, 
do not depend on particular selection of parametrization. 
Moreover, although at the tree level the properties of 
quasiparticles appear to be changed, any correlators in- 
volving the gauge-invariant current retain their values. 
This can be checked by performing the same transforma- 
tion in the presence of additional external gauge field Aq: 
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the reparametrization does not affect the minimal gauge 
couphng with this external field. 

Having established the reparametrization symme- 
try ( p^ of the complete model let us resume the 
transformations of the partition function (0). As usual, 
the field a„p introduced in (^) is defined up to a gauge 
transformation 



(13) 



with arbitrary A„; for simplicity we imply the Coulomb 
gauge A^a„^ — for all our gauge fields. The integration 
over 6n results in 

Z[A] = I ^ — 27rz6„^TO„^ — ?^„^t6„^^ j-^^-j 

or, after regrouping the second term in the exponent, 

where the integer-valued vorticity is defined as 

Mn^l = e^"^ [nin+fi+up - mn+p.p) ■ (16) 



Since the vorticity is locally conserved. 



(17) 



we may exchange the summation over m„^ for a sum 
over Mnfi subject to condition (p7|). This integer- valued 
constraint can in turn be removed by new phases 



(18) 



To finish the transformation of the dual model towards 
the Villain form analogous to the original model (||J|), 
Peskin l40| introduced a convergence factor 



1 = lim exp 



2 -'"n^ 



(19) 



into Eq. ( [l5| ) and transformed the resulting sum 



with the Poisson summation formula 



A/=-oo 



^ gi(/7M-tMV2^ /f^ ^ g-((^-2^m)V2t^ ^21) 



m— — oo 



After rescaling the gauge fields 27ra„p = ^np, 27r6„p = 
_B„p, the dual action can be finally written as 



d6. 



2t 



(22) 



Again, we can integrate away the original gauge field to obtain the final form of the dual action 



Z — lim 



exp 



(23) 



where the appropriate summation in each term of the 
exponent is implied. In accord with our previous conclu- 
sion, the action depends only on the transverse part of 
the combination {JC^^ + TS*)^^. In addition, the dual 
model has a form of (|^,§) with zero lattice temperature 
t = and the dual gauge kernel 

yp2_px^-l p p2 

(2^^ -(^(^ + ^"^ (24) 

where only the transverse part enters because of the 
gauge fixing for the field A. The dual model has the 
same form as the original one and, since the new gauge 
field was introduced as an auxiliary field, this model has 
the same reparametrization freedom (|l2|). Therefore, the 
most general form of the dual model with the lattice tem- 
perature T and the gauge kernel K. satisfies the equation 



f + K:-'^0 + JCo'^^^{JC-'+Ty\ (25) 

or, more symmetrically, 

(r + /C-i)(f + ^-1)==-^. (26) 

An alternative direct derivation of the finite temperature 
dual action, as well as the vortex — monopole mapping 
explicitly relating the scalar sectors of the two models 
is provided in the Appendix ^ This algebraically exact 
mapping proves that dual models of the form (|^,^ related 
by the generalized duality ( p6| ) are equivalent, being just 
two different representations of the same model. 
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To provide an example of the duality in the presence of 
the Chern-Simons coupling, let us consider the simplest 
linear in momenta form of the gauge kernel 



K.'^'ik) = -^e^"'PP„. 
zna 



(27) 



While this definition is very well convergent towards the 
continuum limit and generally looks like a plausible def- 
inition of the Chern-Simons kernel, it is actually non- 
local in coordinate representation. The conventional du- 
ality transformation results in zero lattice temperature 
model (23) with the gauge kernel 



ZTT 



T 



[P^S"" - pfP""] , (28) 



analogous to the frozen superconductor |Q. The two 
terms here are the non-local Chern-Simons term of the 
same form as ( p7|) and the usual lattice Maxwell action. 
The zero lattice temperature implies zero tree level prop- 
agator for the matter field, in analogy with the dual the- 
ory of non-relativistic continuum bosons coupled to 
the purely statistical Chern-Simons gauge field. However, 
already in the RPA approximation, the dynamics of the 
dual gauge field A, associated with non-zero dimension- 
ful charge = (27r)^/T, renders this propagator finite, 
therefore introducing some non-zero lattice temperature 
for the dual matter field. 

Let us try to understand this result using our gener- 
alized duality transformation (p6|). Specifically, we want 
to find a form of the dual model with purely statistical 
gauge field, so that all dynamics would be associated with 
the scalar field. 

To simplify the subsequent algebra, it is convenient to 
introduce special notation for typical matrices appear- 
ing in the expressions. The matrices (5^^ = — PfiPu 
and [I]^^ — e'^'^PPp commute with each other, while 
= —Sj^i^- Algebraically, the symmetric matrix is 
equivalent to unity, while I plays a role of « = exp(z7r/2). 
Since we always assume the transverse gauge, it is pos- 
sible not to write (5* ^ at all and treat I as a commuting 
number. In these notations the dual gauge kernel (|2~ 
can be written as 



Kn = -I- 



,aP TP^ 



27r 



(2n) 



so that the corresponding inverse matrix in the transverse 
gauge is just 




(29) 



Naively, the first term of the expansion may be associ- 
ated with the inverse of the new Chern-Simons kernel, 
while the second constant term plays a role of the dual 
temperature. We notice that while the Chern-Simons 
term is the same as the one in the representation with 
T — 0, the term has been traded for the non-zero 
lattice temperature T — Tjo? in the scalar sector. 

In reality, we cannot just discard the extra terms de- 
noted by ellipsis in Eq. ( p9|) , but we certainly have free- 
dom to choose the value T — Tjo? for the dual tem- 
perature as long as the general duality condition ( [2^ ) 
is satisfied. With this particular choice the exact dual 
gauge kernel becomes 



\aP 



1 



27r 



ia; + (1 + \x) 



(30) 



where we denoted x = TP/2TTa. The exact kernel K, 
differs from the truncated dual kernel 



27r " - 



Ztt 



only by corrections of the higher order AO{P^)A, so that 
in the naive continuum limit the dual gauge field A does 
not have any term or associated dynamics. 

It is tempting to propose that, as long as these cor- 
rections are small, as terms of higher order in momenta 
they should be irrelevant and the scaling properties of 
both the truncated and the exact forms of dual the- 
ory should be equivalent with the corresponding criti- 
cal indices equal. Then, using the previously established 
mapping between theories that are dual as specified by 
Eq. ( p6| ) , one could conclude that the two models with the 
gauge kernels of the form ( |27| ) and the corresponding pa- 
rameters (T, a) and (T/ , — I/a) should be also equiva- 
lent at large scales. The first statement is, however, not 
so obvious because the Chern-Simons theories are intrin- 
sically non-local and the critical dimensions of different 
operators may substantially differ from their classical val- 
ues. We shall prove the legitimacy of this procedure in 



the Section lU 



Along with the dimensionless parameter a we carried 
out the transformation of the classically irrelevant dimen- 
sionful lattice temperature T. Clearly, since the values 
of the lattice momentum P are limited from the above, 
the small values of the ratio T/\a\ = T/\a\ ensure the 
smallness of x, so that the effect of higher order irrel- 
evant terms is negligible and the bare value of the lat- 
tice temperature T should be meaningful. Remarkably, 
in the presence of the Chern-Simons term the duality 
transformation does not lead to reversion of the 

temperature axis characteristic of the duality between 
the superconductor and the superfluid. In our case only 
the Chern-Simons coupling constant is inverted while the 
scalar sector may remain in the strong coupling regime. 

The physical properties of excitations in the general- 
ized x-y model are revealed by their coupling with exter- 
nal fields in 2-|-l dimensions. Luckily, the introduction of 
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the additional external gauge field cAq does not influence 
the exact procedure of the duality transformation; now 
it results in the partition function ( p^ ) up to the substi- 
tution + e^on/^ • Integration over the original 
gauge field A and proper shift of the dual gauge field A 
results in the frozen model (^ ) with additional minimal 
coupling to the external field 

TP 

-Ao, 



1 



and the constant term 



2'Ka 



Ix 



Ao{k) 



(31) 



in the exponent. These two terms remain intact in 
reparametrizcd models with non-zero lattice tempera- 
tures. 

It is clear that the model describes particles with the 
charge e = — e/a, while the momentum-dependent term 
1-f Ix in the denominator of the coupling can be regarded 
as a formfactor indicating the finite size of excitations 
and the presence of the Magnus force. The term ( |3l] ) de- 
pendent only on the external gauge field shows that the 
duality transformation separated the condensate with the 
total charge density — e^A x Ao/27ra = — e^i?o/27ra de- 
termined by the external magnetic field Bq. More specif- 
ically, if the original theory had the average charge den- 
sity p, the total charge density of excitations in the dual 
representation is given by 

One flux quantum of the magnetic field binds the charge 
e/a, or exactly one quasiparticle in the dual represen- 
tation. Since the duality translates the original quasi- 
particles with charge e into vortices with unit vorticity, 
this condensate can be also interpreted as the average 
vorticity 1/a per quasiparticle in the dual model. 

Certainly, the external magnetic field does not at all 
simplify the lattice model: even the one-particle Hof- 
stadter spectrum in the magnetic field is very compli- 
cated. However, if both the magnetic field and the par- 
ticle density are small at the scale of the lattice, only 



the filling factors of electrons v = 2TTp/e^B and quasi- 
particles v = 2'np/eP'B are important. In this case the 
equation (R2h becomes 



i) — a{ah' — 1). 



(33) 



II. 



PERIODICITY IN THE CHERN-SIMONS 
COUPLING 



Non-relativistic continuum models display the period- 
icity in the Chern-Simons coupling because geometrically 
this coupling is just the linking number between the tra- 
jectories of quasiparticles. As long as this linking number 
remains integer, equal increments in the Chcrn-Simons 
coefficient result in the total shift of the phase of the par- 
tition function by integer multiples of 2tt. This is always 
so when quasiparticles avoid each other — for example, if 
they are fermions or if they have strong hard core re- 
pulsion. Although the x-y model is already in the limit 
of point-like strong repulsion as viewed from the corre- 
sponding continuum model, the current lines can inter- 
sect or even join each other and special effort is needed 
to define the Chern-Simons action periodic in its cou- 
pling. In this section we find such definition and use it 
to derive the exact transformation analogous to the flux 
attachment transformation in the continuum, valid for 
an arbitrary form of the lattice Chcrn-Simons action. 

The current of quasiparticles in the x-y model is de- 
fined as the field canonically conjugated with the gra- 
dients of phases A^0„. Since the partition function is 
periodic in these phases, they are cyclic variables and 
the current is integer- valued. As one would expect, the 
current operator is not diagonal in terms of the phases 0„; 
it is more convenient to deal with currents in the special 
current representation. In the Villain approximation the 
explicit form of this representation can be obtained with 
the inverse of the Poisson summation formula The 
integration over phases 0„ ensures the conservation of the 
current Af„^, and the partition function (||,^ rewritten 
in terms of this current becomes 



'T M^^/2 iMniiAn^j, AnfiK-^j^i An' u /2 



(34) 



It is important to emphasize that the transformation to 
the integer current representation ( p^ ) did not change the 
gauge field A in any way, so that the gauge kernel /C^J^/ 
in equations (0) and (p4|) is exactly the same. 



The linking number between lattice current flow lines 
docs not have a natural geometrical meaning at inter- 



section points and therefore cannot be uniquely defined 
to fit all intuitive requirements. There is, however, no 
difficulty in defining the linking number for any two con- 
served fields Af„^ and A^„^ determined on mutually dual 
lattices since they never intersect. Formally, this cou- 
pling can be introduced via the solution aP of equations 
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(35) 



as A/" = J2 ^^npO,^nfi- Since the field Nnp is defined on 
the links of the dual lattice (or at the plaquettes of the 
primary lattice,) the auxiliary field a'^ni^ is defined on the 
links of the primary lattice, and the summation is well 
defined. One can prove that this definition is indeed the 
integer- valued linking number by rewriting the conserved 
current Af„p as a superposition of some number ne of di- 
rected loops L*, i = 1, . . . , carrying unit current each. 
This splits the expression for M into a number of sums 
over independent closed loops and by the Stokes theo- 
rem each of them is exactly equal to the (integer) fiux 
of current N through the surface delimited by the corre- 
sponding loop; obviously the total A/" is the integer- valued 
linking number. 

Now let us define the regularized integer sey-linking 
number of the integer-valued current M in exactly the 
same way but with additional identification Af„^ = Nn^. 
This identifies neighboring parallel links of the origi- 
nal and dual lattices mutually displaced in the direction 
(1,1,1) by half a period. This displacement, uniquely de- 
termining the integer- valued regularized linking number, 
is revealed in Fourier representation, where the auxiliary 
equations (|35|) take the form 



1) 



(k) = e-^'=''/2Mp(k), 



1-e 



— ikj. 



)a°.(k) =0, 



(no summation in p\) or, somewhat more symmetrically, 

e^^Ppy, = ze-'('===+'=«+'=^)/2^fp, P,a°,(k) = 0. (36) 

Solving these equations with respect to auxiliary field a*^, 
we obtain the regularized self-linking number in Fourier 
representation 



M = Tra%(-k)Af^(k) = iTi e'^-'P M^{-k) 



QqPv 



Mp(k), 



where the original exponent exp i /2 fc^ in ( |36[ ) was re- 
placed with Qo = Qo(k) = cosX]^ allowing for the 
symmetry of the total sum. The exponent of the linking 
number Af can be further transformed to form a gauge 
coupling 



exp —inaJV — > 

/ exp^-iM^(-fc)a^(fc) +a(-fc) 

J da ;, 



IF 



(37) 



with new transverse fluctuating field a introduced as the 
Hubbard-Stratonovich field. Clearly, the gauge coupling 
has exactly the form of that in ( 34 ) with the gauge kernel 



Po 



27raQo ' 



(38) 



The momentum expansion of this kernel starts with lin- 
ear antisymmetric term and has no quadratic in momenta 



part; therefore, it corresponds to the Chern-Simons term 
in the long-distance limit. The apparent singularity of 
this gauge kernel at the plane kx + ky + — tt does 
not lead to any divergences but merely suppresses fluc- 
tuations of the field a in the vicinity of this plane. One 
can totally avoid introducing somewhat unpleasant di- 
vision by zero by considering periodic lattices with an 
odd number of lattice nodes in each direction, so that 
at any finite system size the singularity of the kernel is 
never reached. This also eliminates the problem of zero 
functional denominator arising from integration over the 
transverse part of the new gauge field a in (^ . 

Consider the model (^J) with the specific form of the 
gauge kernel (|38|). This partition function is convergent 
at all positive lattice temperatures T even though the 
second term can be zero for certain configurations. One 
can use the Poisson summation formula to rewrite the 
model in the form (|2|j|) depending only on the phases 9. 
By construction, this model is a periodic function of a 
with the period Aa = 2, so that any two models with 
the couplings related by 



a + 2m, T ^ T 



(39) 



are absolutely equivalent to each other. Definitely, this 
periodicity is precisely the flux attachment symmetry as 
it was formulated for non-relativistic systems. It is im- 
portant to mention that none of the previously consid- 
ered definitions [ p8pl| -|3^ of the lattice Chern-Simons 
term provide for this property; this is why we needed 
to construct yet another form of the Chern-Simons cou- 
pling. The flux attachment transformation changes the 
properties of the gauge field, but our main interest is to 
understand the effect of the fractional statistics on the 
dynamics of the scalar field; we use the auxiliary Chern- 
Simons field only to get rid of the non-local interaction. 

Even though one can construct other forms of the 
gauge coupling leading to the same symmetry of the 
partition function, these forms are very special; generic 
gauge couplings do not reveal the exact fiux attachment 
symmetry. Moreover, it is easy to check that this prop- 
erty is not even preserved by the duality transforma- 
tion (^6|). Therefore, we need to define the flux attach- 
ment transformation to display the same periodicity ( |39| ) 
at least in the continuum limit for some broad class of 
possible lattice Chern-Simons terms. 

In order to do this, let us introduce the trivial phase 
2T:mJ\f proportional to the integer-valued self-linking 
number A/" in addition to the already present in Eq. ( |3^ ) 
gauge coupling with the gauge kernel /C. Clearly, this 
trivial phase does not change the partition function and 
yet, after rewriting the introduced term as the addi- 
tional gauge coupling ( ^7[ ) and integrating away one of 
the gauge fields, we obtain the same model but with the 
gauge kernel 



P 



(40) 
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The resulting model is rigorously equivalent to the orig- 
inal model with the gauge kernel /C and it is convergent 
at any non-zero lattice temperature T. 

The long-range properties of the flux-attached gauge 
kernel (^) depend on the properties of the original ker- 
nel /C. If /C has the Chern-Simons form, i.e. its momen- 
tum expansion starts with the linear in momenta term 
l^fiiy ~ the expansion of the new gauge ker- 

nel JC has the same form with the coefficient a' — a-f 2m. 
In the particularly simple case when the matrix struc- 
ture of the gauge kernel /C^^ is the combination of i5* ^ 

and [I]^i/, we may further simplify the flux attachment 
transformation by defining the formfactor Q through 
/C = IP(5~^/27ra. Then the appropriately defined form- 
factor Q' for the new kernel /C' is determined simply as 
the linear combination 



Q' = 



(41) 



Earlier we adjusted the lattice temperature to obtain the 
purely statistical Chern-Simons field without any dynam- 
ics in the naive continuum limit, the appropriate lattice 
gauge kernel having no quadratic terms in small momen- 
tum expansion. This is equivalent to the expansion of 
the formfactor Q(fc) = 1-1- O(fc^), fc ^ 0; obviously this 
property is preserved by the equation (^l|). 

Before discussing the implications of this flux attach- 
ment transformation on the universality of the phase 
transitions in Chern-Simons models, let us reflect on the 
fact that the derived expression ( ^0[ ) is an exact sym- 
metry of any lattice model with the integer-valued cur- 
rent, moreover, it is valid not only for Villain models, 
but also for usual x-y models minimally coupled with 
with some fluctuating gauge field. If the original model 
has no long-range interaction between the quasiparticles, 
it is in the universality class of the x-y model, while the 
transformed model has a Chern-Simons term with the 
gauge kernel ( |38| ) and the coupling a = 2m. The duality 
transformation results in the model of the same form with 
the lattice temperature T — T jArn? and Chern-Simons 
coupling a = — l/2m, the small momentum expansion 
being valid as long as TP <C 47rm. The continuum limit 
of this model is that of scalar particles coupled with the 
even-denominator Chern-Simons field; yet this theory is 
rigorously equivalent to the original x-y model! 

Completely different situation arises if we try to add 
such additional Chern-Simons interaction to lattice su- 
perconductor described by the same model but with 
quadratic in momenta gauge kernel /C. Qualitatively, the 
interaction introduced with the CS coupling is working 
independently of the original one. In coordinate represen- 
tation this additional interaction between current lines 
decays as l/r at large distances; in order to have any ef- 
fect it should not be screened by the original interaction. 
However, the current lines in normal state already have 
the magnetic interaction logarithmic at large distances, 
and the additional coupling cannot have any effect on the 



long-distance physics: it will not be visible at all at some 
finite scale. 

Although this statement appears to be rather evi- 
dent, in some cases the form of interaction at large dis- 
tances may not be obvious. For example, Schultka and 
Manousakis considered the 3D x-y model with the 
Chern-Simons coupling introduced for vortices instead 
of the real charges. Numerically, the phase transition in 
this model is in the x-y universality class for all values 
of the Chern-Simons coupling, and this was interpreted 
as the explicit proof of universality of phase transitions 
in the Chern-Simons theories. We do not agree with this 
interpretation: clearly, the model dual to the x-y model 
is the frozen lattice superconductor and the CS coupling 
between the vortices of the original model is just the ex- 
tra coupling between the currents of the dual theory. But 
the currents in the normal state of superconductor al- 
ready have the long-range magnetic interaction, the ad- 
ditional coupling is screened, and effectively the large dis- 
tance behavior is merely that of the original x-y model. 
This statement can be obtained more rigorously from the 
small-momentum expansion of the combined gauge ker- 
nel derived in the Appendix pi 



III. THE LAW OF CORRESPONDING STATES 



For the Villain x-y model with fractional statis- 
tics a introduced by coupling to an auxiliary gauge field 
with the kernel obeying K.^"' = e^^'Pkpl^-Ka 0(k^) at 
small momenta, we found that the exact duality trans- 
formation (|26l T = T jo?) and flux attachment transfor- 



mation (Ml T' — T) have very natural and simple form 

a = — 1/a, T ~ Tjc? 
a' = a + 2m, T' ^ T. 

in terms of the Chern-Simons coefficient a. Although 
a fully describes the gauge coupling only in the limit of 
small momenta, we would like to check whether these 
relationships are meaningful by themselves at least for 
some values of the lattice temperature T and the Chern- 
Simons coupling a. 

In the continuum, the duality and fiux attachment 
transformations can be used to "unwind" the Chern- 
Simons coupling with the fractional parameter — 
Po/qo. Shifting ag by an appropriate even integer num- 
ber, we can always satisfy inequalities —1 < o'q ~ 
ao -I- 2mo < 1- If aj, = 1, or af, = 0, we can- 
not proceed any further with our transformations and 
should stop. Otherwise the duality transformation yields 
ot\ = — l/cto = Pi/qi with new integer denominator 
< gi < go- Repeating such transformations, in a fi- 
nite number of steps n < qo we arrive at either a'^ = 1 
or Qf^ = with the final lattice temperature r„ — q^T^ 
given simply by the denominator go of the original Chern- 
Simons coupling. Since both the shift and the duality 
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transformation preserve the parity of the sum of numer- 
ator and denominator of a, it is clear that any fraction 
with this sum odd (i.e. either numerator or denominator 
is an even number), will lead to a„ = 0, and the original 
fraction with both p and q odd will result in a„ = 1. 

In the former case the final theory is just the pure x- 
y model without any additional coupling. Up to finite 
renormalization correction, the phase transition temper- 
ature in this model is just that in the usual Villain x-y 
model [|l|, r„ - TvxY = 3.03. We immediately obtain 
values 



To a = — 



qo 



(42) 



of the phase transition temperatures in the Villain Chern- 
Simons models with either numerator or denominator 
given by even numbers, as plotted in Fig. |l| with solid 
lines. Similarly, once the phase transition temperature in 
a single x-y model with "fermionic" Chern-Simons cou- 
pling a = 1 is known, one can use the same relationship 
to obtain the approximate values of the phase transition 
temperatures for all "fermionic" fractions with both nu- 
merator and denominator odd. 

Let us rectify the outlined procedure by using the ex- 
act duality ( p6| ) and the flux attachment (^0|) transfor- 
mations at every step, starting with the local theory with 
the gauge kernel (|3^) and the initial parameters ao and 
Tq. The partition function with this gauge kernel is pe- 
riodic with respect to the coupling ao, and we can safely 
assume that < |ao| < 1. The dual model has the lattice 
temperature Ti = Tq/uq > Tq and the gauge kernel 



IP 



1 



2iTai 



Ixq + {Qo + 



(43) 



where xq = Pro/27rao, \xo\ < To/n\ao\ and the the new 
coupling ai = — 1/ao. It is convenient to rewrite the new 
gauge kernel as 



IP 



2naiQi ' 

by introducing the next level formfactor 

1 



Qi Ia;o 



Qo + Ixq 



(44) 



(45) 



Similarly introduced formfactor of the CS coupling re- 
sulting from the flux attachment transformation (^0|) can 
be expressed merely as 



aiQi + 2miQo 
ai + 2mi 



(46) 



where ai + 2mi = a'l is the new effective gauge coupling; 
obviously, both transformations preserve the small mo- 
mentum expansion of the formfactor Q{k) = 1 + 0{k^). 



At the n-th step a„ becomes an integer and the univer- 
sality class of the resulting theory is determined by the 
parity of a„. 

If this number is an even integer, the final shift simi- 
lar to ( ^ ) results in the gauge kernel singular at small 
momenta, and the long-distance gauge interaction is sup- 
pressed. The remaining short-range interaction between 
the currents perturbs the x-y critical point, and one can 
prove that it is irrelevant in this point by simple power 
counting. As formally irrelevant, this additional inter- 
action may result either in some finite correction to the 
transition temperature, or it has to change the symmetry 
of the phase transition completely. In the former case, 
since the models at every level of hierarchy are exactly 
equivalent to the original Chern-Simons model with the 
fractional coupling ao, we conclude that the phase tran- 
sition in the original model is in the same x-y universality 
class. 

Similar arguments apply when the final a„ is an odd 
number corresponding to the Fermi statistics of the 
quasiparticles. Although we do not know the universal- 
ity class of the phase transition for Chern-Simons bosons 
with odd-integer-valued coupling, we can claim that it 
should be the same for all original models in this class 
unless the irrelevant terms drive it away from the critical 
point. 

So far we concentrated our attention on more pleas- 
ant possibility that the irrelevant terms are too weak to 
change the symmetry of the critical point and the phase 
transition does not change its universality class or be- 
come the first order. Let us try to analyze the irrelevant 
terms more carefully to see what fractions are likely to 
follow the universality scenario. 

We saw that the bare lattice temperatures of the tran- 
sitions decrease rapidly with the hierarchy level, there- 
fore the temperature-dependent convergence condition 
Xk = PTk/2'Kak <SC 1 is not very limiting. There is, how- 
ever, the momentum contribution from the formfactor 
Qo = cos^^ kfj/2 of the local Chern-Simons kernel (|3^). 
To investigate the extreme possible effect of this contribu- 
tion, let us write the sequence of hierarchical formfactors 
at zero lattice temperature for the fraction 



ao = 



(47) 



2mi 



2m2 



starting with the formfactor Qo- The first duality trans- 
formation in the unwinding procedure results in Qi — 
1/Qo, Oil — — 1/ao, and the flux attachment transforma- 
tion (^) leads to the formfactor 



1 

Q'o 



1 — 2toi 2m2 



1 



2m3 + 



V 



where 



V' = l-Qi = 



2 =sin2^fc^/2 



(48) 



(49) 
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vanishes at the origin. Clearly, the coefheient in front 
of the quadratic in momenta part increases with the in- 
creased levels of hierarchy, and eventually it may become 
the main driving term of the phase transition in the sys- 
tem. 

Making yet another duality transformation we obtain 
Q2 = l/Q'i! the new zero-temperature kernel has zero 
surface at finite momenta determined by the equation 
-P^ = l/2mi(2TO2 -I- (l/2rn3 + ...))< l/2mi. Although 
the divergence does not occur at any finite lattice tem- 
perature, the system develops a soft mode much closer 
to the origin then the location of the original soft plane 
in the middle of the Brilluen zone — I. The exact 
location and the orientation of this soft surface depends 
on the details of the selected regularization procedure (in 
our case it was determined by the chosen form of the local 
kernel (^8|)), but its existence is probably unavoidable as 
long as we want to define the integer- valued linking num- 
bers. The fluctuations in the vicinity of this mode will 
grow with the hierarchy level and, again, are capable of 
destroying the second order phase transition. 

We see that the "unwinding" procedure creates ad- 
ditional instabilities of the strength increasing with the 
level of hierarchy and, generally, with the denominator 
of the statistical coupling. Therefore, the phase tran- 
sitions in the system of particles with fractional statis- 
tics a — p/q are expected to be universal only for small 
enough denominators q. 

One could argue that instead of starting with the 
Chern-Simons bosons with fractional statistics, we could 
have followed the usual direction of the hierarchy se- 
quence and started with the x-y model or Chern-Simons 
bosons with odd "fermionic" statistical coupling. In this 
case the sequence of exact transformations also results in 
non-local Chern-Simons model with fractional coupling, 
and the truncated local model with the same Chern- 
Simons coefficient differs from the non-local one by clas- 
sically irrelevant terms. Again, the two models should be 
in the same universality class as long as the phase tran- 
sitions remain of the second order. This procedure has 
the advantage that one can construct the sequence of lat- 
tice models without any apparent divergences. Indeed, 
starting with the ao — 2mo, the duality ( [l5|) and the fiux 
attachment (|46|) transformations result in the formfactor 



1 



( 



V 



(50) 



2mi 



2mn 



the coefficient in front of here is less then one and 
the original convergence condition \k\ <C 1 is preserved. 
Unfortunately, this argument also fails to prove the uni- 
versality at high hierarchy levels: the perturbation is not 
small compared to the relevant scale given by the van- 
ishing lattice temperature (p2[). 



CONCLUSIONS 



In this paper we studied a relativistic version of the 
Chcrn-Simons-Landau-Ginzburg theory of bosons in the 
limit of strong coupling. We used the lattice representa- 
tion of this model in terms of Villain x-y model minimally 
coupled with the Chern-Simons gauge field to access the 
strong coupling limit without the perturbation theory. 
This model has duality and flux attachment symmetries 
in the long wave length limit, i.e. these symmetries are 
accurate up to irrelevant cubic and higher order deriva- 
tive terms. There is no single lattice model that can obey 
both symmetries exactly, but we construct algebraically 
exact non-local duality and flux attachment transforma- 
tions corresponding to these symmetries in the contin- 
uum limit. These non-local transformations were used 
to to show that there arc only two universality classes in 
this model, one corresponding to the pure x-y transition, 
and another corresponding to "fermionic" x-y transition, 
or the transition in the x-y model with the Chern-Simons 
coefficient a = 1. The value of such an investigation is 
to establish a theoretical model in which the universality 
in a CSLG theory can be proven beyond the bounds of 
perturbation theory. 

Although it may be possible to construct an experi- 
mental system that would be described by such relativis- 
tically invariant model, this model is not in the same 
universality class as the experimental quantum Hall sys- 
tems, where the relativistic symmetry is broken by non- 
zero charge density, external magnetic field and the dis- 
order potential. Theoretically, these effects can be incor- 
porated in the lattice model as additional external gauge 
fields without breaking the exactness of the performed 
hierarchy transformations, but the extra flelds lower the 
symmetry of the problem. Now every phase is character- 
ized by the non-trivially transforming filling factor v in 
addition to the original Chern-Simons coefficient a, and 
the exact mapping between different quantum Hall states 
is absent. 

This is not surprising, since the mapping preserves 
the information about the exact configuration of dis- 
order. We believe, that the disorder averaging with 
proper identification of the relevant degrees of freedom 
will increase the symmetry and reveal the universality of 
phase transitions in this model. Qualitatively, we saw 
that the Chern-Simons interaction has no effect if it is 
screened by some other independent long-range force. 
Two-dimensional time-independent random scalar poten- 
tial may be treated as the interaction of infinite range 
in time direction; apparently this interaction is strong 
enough to suppress or modify the effect of the statistical 
coupling on localization phase transitions. 

Our model also sheds some light on the important 
question of existence of high-order hierarchy states in the 
quantum Hall effect. Since the duality and periodicity 
transformations can simultaneously become symmetries 
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of the theory only in the long-distance limit, the accumu- 
lation of extra terms, irrelevant near the critical point, 
will eventually render the phase transition of the first 
order or change its symmetry. 

Further research on this class of models should be 
concentrated on understanding the models with several 
species of the scalar fields. Clear understanding of the 
model with finite number of components is necessary to 
apply the replica trick, which is the only possible way to 
perform disorder averaging in the presence of interaction. 

Another perspective direction is to understand the hi- 
erarchical picture within the rigorous approach by Chen 
and Itoi . Starting from the representations of the 

2-1-1 dimensional Lorenz group, they established the ex- 
act relationship of the spin of particles with the coefficient 
of their Chern-Simons coupling, namely that all theories 
of particles with spin s and Chern-Simons coupling a de- 
scribe the same irreducible representation of the Lorenz 
group with the fractional spin s + a/2. Although nomi- 
nally the action of the particles with the spin s has 2s H- 1 
components, in 2 + 1 dimensions the equations of motion 
for massive particles leave only one independent compo- 
nent, and the long-range properties of all particles with 
the same a but the spins s differing by an integer should 
be similar. This can be considered as the rigorous defini- 
tion of the flux attachment transformation for hard-core 
relativistic particles with Chern-Simons interaction. 

While finishing this work we learned that a similar 
model is studied by Fradkin and Kivelson |5^. They 
also use duality and periodicity transformations to ar- 
gue that phase transitions in three-dimensional lattice- 
regularized models are universal, although they are less 
specific about the form of the Chern-Simons coupling at 
short distances. In addition to the behavior considered in 
this paper, Fradkin and Kivelson consider a model with 
additional dimensionless non-local interaction between 
the integer- valued currents, corresponding to the dissipa- 
tive conductivity in the quantum Hall effect. As in Refs. 
1^ , ^ , this leads to the phase diagram with multiple 
self-dual fixed points. Perturbatively, such interaction 
of finite strength can be generated near the phase tran- 
sition point ||2^ , and the accumulation of higher-order 
irrelevant terms considered in the present work can in 
principle imply the change of the symmetry of the sec- 
ond order phase transition instead of the first order phase 
transition. 
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APPENDIX A: THE ALTERNATIVE 
DERIVATION OF THE FINITE TEMPERATURE 
DUAL MODEL 

Instead of using the previously established symmetry 
property ([l^) , we could have arrived directly at the dual 
action with non-zero lattice temperature by equivalent 
transformations of the gauge field in the exponent. In- 
deed, averaging the equation (|l^) in the original gauge 
field Anfj. produces the kinetic term (^J) for the dual 
gauge field a, and one can further rewrite the exponent 
in terms of new auxiliary transverse field b as 



iaM 



1 



1, 



aK-^a iaM — iab &/Cg 



Now one can add and subtract the term with 
the arbitrary constant T, so that the introduction of 
yet another transverse field c to dispatch off the part 

b (JCq^ — 6/2 of the obtained gauge kernel results in 
the exponent 



iaAI — iab + ibc 



T. 



i.(*o-' 



This expression is linear in the field a and an easy in- 
tegration yields the constraint b = M; the subsequent 
integration in b trivially yields 



iMc- -M^ - 
- 2 



T 



With the transverse integer- valued current M this is pre- 
cisely the exponent of expression (|2^) with the gauge ker- 
nel K. = (j^Q^ ~ and the new lattice temperature 

T > replacing the artificially introduced infinitesimal 
variable t. Now the summation formula (21) and subse- 
quent rescaling of the new gauge field lead directly to the 
dual model of the form (|^,^ with new finite temperature 
T and the gauge kernel JC satisfying (^). 

The derived mapping is a generalization of the well- 
known duality in three dimensions [^ . It can be un- 
derstood as the duality between vortices and monopoles 
in three spatial dimensions, or the charge — vortex dual- 
ity in 2 -|- 1 dimensional systems. Indeed, the average in 
the presence of an arbitrary number of vortex-antivortex 
pairs introduced by the integer vortex strength L„ 



exp i 



exp I 



I^l_.a) (A1) 



with zero total vorticity J2n ~ ^ results in the con- 
straint 



^^2L^9n - l6rA,bn, ^ g ( ^^^^^^ _ ) (^2) 



instead of (||). Evidently, the integers L„ serve as the 
sources for the field b, or the monopoles with appropriate 
quantized charges. The values of these charges are not 
affected by any single- valued transformation of the vector 
potential, therefore the statement of equivalence gener- 
ally holds for any two models related by Eq. (|2^). Simi- 
larly, the dual transformation of the original model (||J|) 
in the presence of monopoles introduced by appropriate 
multivalued external vector potential results in vortices 
of the form (Al) located exactly in the original positions 
of the monopoles. 

Unfortunately, any averages involving the gauge fields, 
particularly the averages of gauge-invariant currents, do 
not have simple dual representation since the ffuctuating 
gauge fields change upon reparametrization. Therefore, 
only scalar sectors of arbitrary models of the form (|^J|) 
related by the generalized duality transformation (^6|) are 
equivalent to each other. 



APPENDIX B: CS COUPLING BETWEEN THE 
VORTICES OF X-Y MODEL: DETAILED 
ANALYSIS 

Here we analyze the model numerically considered by 
Schultka and Manousakis [p2|: the pure x-y model with 
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the additional Chern-Simons interaction between the vor- 
tices, introduced to study the universality of transitions 
in a system of particles with fractional statistics. The 
original model is defined simultaneously in terms of orig- 
inal phases 9n and the vorticity M„^. Working in the 
Villain approximation, let us perform the exact duality 
transformation on the scalar sector and express the model 
entirely in terms of dual variables. The duality transfor- 
mation results in the model of the form ( pO[ ) with the ad- 
ditional coupling of the vorticity M to the second gauge 
field A with the Chern-Simons kinetic term. Shifting the 
gauge field 

and integrating away the field a we obtain the gauge ker- 
nel of the combined interaction in the form 



IP 
2TTaQ 



1 



i^"^ 1 - IT&QP/2TT J 
0{P^), aTPQ < 27r, 



(Bl) 



47r2 

instead oi T P"^ / Air"^ present without any gauge coupling; 
the statistical coupling of the vortices does not change 
the long-range properties of the model at all! 

This is precisely the result of numerical study [ p2[ ; 
from our analysis it is clear that the additional Chern- 
Simons interaction between the vortices gets completely 



screened by the Coulomb interaction already present in 
the dual model. Indeed, the original vortices of the x-y 
model map to charges in superconductor, which is the 
correct dual model. The charges have the Coulomb in- 
teraction decaying like 1/r at large distances, while the 
additional CS interaction falls off like l/r^; there always 
exists some range above which the usual Coulomb force 
dominates. It is because of the existence of this finite 
range the asymptotic form is correct only in the restricted 



range of momenta as specified in (Bl) 



This phenomenon can also be understood in terms of 
the original x-y model: the vortices in superfluid are 
not really local objects because of their long-range phase 
structure, therefore the additional Chern-Simons inter- 
action does not do anything at large enough distances. 

It is important to emphasize that the previous analysis, 
including the equation (Bl), holds for x-y model coupled 
to a Chern-Simons gauge field via the dual current only. 
If the long-range structure around the vortices is already 
broken by some mechanism, the addition of the Chern- 
Simons interaction between vortices does change the pa- 
rameters or even the qualitative behavior of the model. 
In this case the long-range phase structure of the vortices 
is destroyed (in the dual representation — the charges are 
screened) and the additional statistical coupling is the 
main interaction at large distances. 



3 



o 
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-2-1012 
Chern-Simons coupling a 

FIG. 1. The phase diagram of the Chern-Simons Villain x-y model. Vertical lines represent ordered phases of different 
symmetry, corresponding to different levels of hierarchy. The phases with the sum of numerator and denominator odd can be 
mapped to usual Villain x-y model, they are shown with solid lines. The phases with "fermionic" Chern-Simons coupling are 
only schematically shown with dotted lines since we know only the ratios of corresponding transition temperatures. Unlike 
the case of usual superconductor, duality transformation does not invert the direction of the lattice temperature axis, and the 
high-temperature phases are always disordered. 



